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PROPAGATION  OF  CHARGED  PARTICLE  BEAMS 
IN  THE  ATMOSPHERE 

Introduction 

This  paper  presents  a  brief  review  of  intense  relativistic  electron 
beam  propagation  in  "dense"  neutral  gas.  To  specify  the  meaning  of 
"dense",  we  shall  assume  that  the  beam  energy  is  deposited  locally  and 
promptly  in  the  gas,  resulting  in  the  formation  of  a  highly  collisional 
plasma  with  electron  density  ng  low  compared  to  the  gas  atom  density,  but 
high  compared  to  the  beam  density.  These  assumptions  are  usually 
reasonable  for  gas  densities  above  a  few  tens  of  torr,  on  up  to  atmospheric 
density.  This  work  was  originally  prepared  for  the  annual  Particle 
Accelerator  Conference  in  Washington,  DC,  March,  1987,  and  was  constrained 
in  length  by  the  requirements  of  that  conference.  It  is  by  no  means  a 
complete  review  of  the  subject,  but  I  have  attempted  to  present  at  least  a 
coherent  overview.  The  work  reviewed  is  the  product  of  many  investigators; 
little  or  no  new  work  is  presented  here. 

A  charged  particle  beam  injected  into  a  neutral  gas  begins  immediately 
to  ionize  the  gas.  For  example,  in  air  at  standard  atmospheric  density 
each  beam  electron  collisionally  ionizes  about  100  thermal  electrons  per 
cm.  Thus,  the  "plasma"  electron  density  ng  exceeds  the  beam  density  n^ 
almost  immediately.  Over  a  time  scale  1/4 no,  a  very  strong  radial  space 
charge  field  expels  a  small  fraction  of  the  plasma  electrons,  thereby 
setting  up  a  charge-neutral  region  within  the  beam  and  out  to  a  large 
radius  b  =  c/4na.  (Since  the  plasma  is  typically  weakly  ionized  and 
collisional,  it  can  be  characterized  by  a  local  scalar  conductivity  a.) 
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Thereafter,  the  only  radial  forces  on  the  beam  are  magnetic:  the  beam  is 
pinched  by  its  self-force,  but  this  may  be  partially  neutralized  by  reverse 
currents  induced  in  the  plasma.  The  beam  thus  propagates  in  a  self-pinched 
equilibrium,  with  the  magnetic  pinch  balancing  the  beam's  transverse 

pressure,  as  well  as  any  centrifugal  force  if  the  beam  is  rotating.  The 
equilibrium  radius  is  proportional  to  the  emittance  and  inversely 

proportional  to  the  square  root  of  the  net  (beam  plus  plasma)  current  In- 
{If  the  radial  profile  of  the  plasma  current  J^fr)  differs  from  that  of 
beam  current  J^r)  it  is  necessary  to  define  an  appropriate  radially- 
averaged  "effective  current"  that  controls  the  pinch  strength.'*'  We  shall 
neglect  such  subtleties  here.}  Since  I  generally  increases  as  one  moves 
backward  in  the  beam,  the  equilibrium  has  a  trumpet-like  appearance,  with 
the  radius  steadily  decreasing  (Fig.  la).  Indeed,  the  very  front  of  the 
beam,  where  charge  neutrality  has  not  been  established,  is  unpinched,  has  a 
large  radius  and  constantly  erodes  due  to  both  radial  expansion  and  energy 
loss . 

Single-Particle  Energy  Loss 

Highly  relativistic  electrons  lose  energy,  due  to  ionizing  collisions 

2  3 

with  gas  atoms,  at  a  rate  given  by  Bethe's  formula,  ’ 

dE/dz  =  -(2rmZe4/mc2)ln(y3m2c4/2tf2<«>2),  (1) 

where  n  is  the  gas  atom  number  density,  Z  the  atomic  number,  and  h<w>  a 
characteristic  bound  electron  energy.  For  electrons  with  E  >  1  MeV  in 
full-density  air,  this  stopping  power  is  200  to  300  keV/m. 

For  high-energy  electrons  (e.g.,  >  100  MeV  in  air),  energy  loss  due  to 
bremsstrahlung  emission  dominates.  This  energy  loss  process  proceeds 
exponentially,  with  the  mean  energy  <E>  decreasing  as 

d<E>/dz  =  -<E>/>  ,  (2) 


iA] 
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with  the  radiation  length  given  hy 

1/Xr  =  4nZ(Z+l)(e2/tfc)(e2/mc2)2lnA,  (3) 

1/3 

with  A  s  192/Z  if  r  >  100,  or  A  s  y  if  y  <  100.  In  standard-density 
air,  X^  =  300m.  Bremsstrahlung  energy  loss  is  statistical  in  nature, 
leading  to  a  large  energy  spread  ("straggling").  This  has  important 

4 

implications  for  beam  stability  and  range. 

Other  single-particle  energy  loss  mechanisms,  e.g.,  synchrotron 
radiation,  are  relatively  unimportant  for  beams  propagating  in  air. 


Radial  Expansion 

As  a  result  of  multiple  small-angle  scattering  off  nuclei,  an 
unpinched  electron  beam  in  standard-density  air  expands  as  the  3/2  power  of 
distance.  This  expansion  is  rapid  as  compared  to  energy  loss.  A  pinched 
beam  also  expands  as  its  emittance  increases  by  scattering,  but  in  a 
different  and  much  slower  way.  Energy  loss  also  has  an  effect  on  the  beam 
radius  a. 

Let  us  consider  a  beam  which  is  subject  to  the  energy  loss  mechanisms 
discussed  above,  as  well  as  to  multiple  small-angle  scattering  at  an 
angular  rate2 

S  s  d<02>/dz  =  16rm(Ze2/Ymc2)2ln(210/Z1/3).  (4) 

Because  neither  electron-electron  collisions  nor  bremsstrahlung  emission 

result  in  angular  scattering  of  a  high-energy  electron  (to  within  order 

1/y,  assumed  negligible),  it  is  easily  seen  that  these  mechanisms  leave  the 

2 

unnormalized  emittance  e  =  a-IO  >  invariant.  If  we  assume  that  scattering 

is  slow  compared  to  a  betatron  oscillation  wavelength  Xa  of  the  beam 

p 

2 

electrons  in  the  pinch  potential,  i.e.,  that  XQS/<0  >  <<  1,  it  can  be 

P 

shown3  that  the  beam  radial  profile  assumes  a  self-similar  Bennett  profile, 

2  2  2 
and  that  £  increases  at  a  rate  a  S.  In  equilibrium,  <8  >  is  fixed  by 

Bennett  condition, 


(5) 
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<e  >  =  i  /i., 

n  A’ 


where  1^  s  17(3y  kA  is  the  Alfven-Lawson  current.  Thus,  the  slow  changes  in 
y  (energy  loss)  and  e  (scattering)  lead  to  a  steady  adiabatic  change  in  a, 

-1/2  1/2 

d  ln( y  i/zIni/za)/dz  =  1/L^,  (6a) 


where 


Lm  s  0ic/2ite2)(I  /I . )y2X  . 
N  n  A  r 


For  example,  this  gives  =  (InE/7xlO  watts)Xf  in  standard  density  air. 
Equations  (6), (7)  are  known  as  the  Nordsieck  equation,  derived  in  its  basic 
form  (but  never  published)  by  A.  Nordsieck  in  the  early  1960s.  Later 


contributions  were  made  by  Lee. 


Fawley  (unpublished  recent  work)  was 


the  first  to  derive  the  correct  energy  dependence  in  the  LHS  of  (6a). 

Energy  loss  that  is  due  to  the  effect  of  E electric  fields  (ohmic 

loss,  to  be  discussed  below)  conserves  the  transverse  momentum  p^  =  ymv^  of 

2 

each  beam  electron,  and  thus,  increases  the  transverse  pressure  <yn^mv^  >. 
Ohmic  energy  loss  thus  causes  beam  expansion,  whereas  it  is  evident  in 
Eq.  (6a)  that  single-particle  energy  loss  leads  to  beam  contraction. 
Mathematically,  it  can  be  shown  that  ohmic  loss  leaves  the  normalized 
emittance  ys  constant,  and  that  when  it  predominates  the  Nordsieck  equation 
takes  the  form 

1/2  1/2 

d  ln(y  In  a)/dz  =  1/Lj,.  (6b) 

A  few  electrons  undergo  larger-angle  single  scatterings,  which  lead  to 
"Moliere  scattering".2  These  electrons  escape  to  large  radius,  and  it  is 
simplest  to  regard  them  as  lost.  This  leads  to  a  slow  decrease  in  beam 
current  1^,  but  to  a  reduced  rate  of  emittance  growth  for  the  remaining 

O 

beam.  The  net  effect  is  to  increase  Eq .  (7)  for  L^  by  20%  to  40%,  which 

9 

is  in  good  agreement  with  experiment. 
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A  beam  can  only  be  considered  to  be  pinched,  and  the  Nordsieck 
equation  only  applies,  if  the  Nordsieck  expansion  rate  (basically  an 


exponential  process  with  e-folding  range  L^)  is  slower  than  the  diffusive 
expansion  rate  for  an  unpinched  beam.  This  sets  a  minimum  condition  for 
propagation  in  the  pinched  mode  in  a  dense  medium.  It  is  also  to  be  noted 
that  in  air  Nordsieck  expansion  occurs  more  rapidly  than  single-particle 
energy  loss  if  I  <  80kA  and  E  <  lOOMeV,  or  if  I^E  <  7x10^  watts  and  E  > 
lOOMeV.  It  is  thus  evident  that  high  current  and/or  energy  are  necessary 
for  effective  self-pinched  propagation  in  dense  media. 

Energy  loss  and  radial  expansion  represent  fundamental  limitations  on 
range,  which  can  be  alleviated  only  by  propagating  in  a  reduced-density 
channel  ("hole-boring").  This  can  be  accomplished  by  sacrificing  the  front 
of  a  particle  beam  to  heat  the  air  and  induce  radial  expansion,  which 
increases  the  range  of  the  beam  tail,  or  by  using  a  series  of  beam  pulses 
to  produce  the  same  effect,  or  by  using  some  other  means  to  heat  and 
prepare  a  reduced-density  channel. 

Plasma  Return  Current,  Ohmic  Loss,  and  Nose  Erosion 

Lenz's  law  suggests  that  a  CPB  should  induce  an  axial  electric  field 
Ez  which  opposes  the  propagation  of  the  beam  and  drives  a  reverse  current 
in  the  conducting  medium.  Ez  extracts  energy  from  the  beam  and  eventually 
dissipates  it  in  the  plasma  through  resistive  decay  of  the  return  current. 
For  a  highly  relativistic  beam  with  v^  =  c,  does  not  arise  at  the  very 
front  of  the  beam,  where  the  gas  is  nonconducting.  The  fields  there  are 
purely  transverse  electrostatic/magnetostatic.  In  effect,  the  beam  serves 
as  a  guide  for  an  electromagnetic  wave  (E  ,B_)  in  the  vacuum.  As  a 
increases,  many  things  happen  in  rapid  succession  at  the  point  where 
C  =  ct-z  =  c/4rta.  (C  is  the  distance  behind  the  beam  head,  a  very  useful 
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coordinate  for  many  purposes.)  The  beam  charge  is  neutralized  and  the 
self-pinch  is  established.  (Hence,  this  region  of  the  beam  is  called  the 
"pinch  point".)  Maxwell's  equations  reduce  to  Ampere's  law  out  to  a  large 


radius  b  wherein  space  charge  neutrality  prevails.  The  E  field  "turns 

around"  and  becomes  an  field  which  is  governed  by  the  inductive  term  in 

Ampere's  law.  This  leads  to  a  very  large  spike  in  the  field  Ez(z),  which 

can  reach  many  MV/m.  As  a  continues  to  increase  rapidly,  the  "monopole" 

2 

magnetic  decay  length  ctq  s  (2noa  /c)ln(b/a)  becomes  much  larger  than  the 

beam  radius  a  and  the  net  current  is  frozen  in.  Thereafter  E  takes  the 

z 

value  necessary  essentially  to  maintain  the  value  of  I  established  at  the 
pinch  point.  Since  a  is  rising  rapidly,  EZ(C)  decreases  rapidly,  making 
the  Ez  spike  very  narrow.  Thus,  ohmic  energy  loss  extracts  energy 
primarily  from  beam  electrons  near  the  pinch  point.  Furthermore,  ohmic 
loss  results  directly  in  radial  expansion,  as  we  have  seen.  Thus,  it  is 
appropriate  to  regard  ohmic  loss  as  primarily  a  mechanism  for  erosion  of 
the  beam  front.  If  one  assumes,  for  convenience,  that  the  beam  current 
Ib(C)  and  voltage  V  are  constant,  then  the  ohmic  energy  loss  rate  is 
equivalent  to  erosion  of  the  beam  front  at  a  rate^’^ 

dC/dz  =  (In/IA)ln(b2/a2).  (8) 

Ve  note  that  nose  erosion  is  additionally  driven  by  scattering,  since 
Nordsieck  expansion  is  fastest  at  the  pinch  point,  where  the  pinch  force  is 
weak.^  Scattering-driven  erosion  is  not  included  in  (8). 

Calculations^  also  show  that  if  the  ionization  of  the  gas  is  due 
entirely  to  beam  collisions  with  gas  atoms,  then  typically  I  settles  down 
to  a  slowly  varying  value 

In  =  Ib/( 1+X)  (9) 


2 

shortly  behind  the  pinch  point.  Here  X  s  dt^/dt  s  d(rioa  /2c)/d£  is  a 

normalized  measure  of  the  beam  current.  The  fractional  current 

neutralization  f  =  -I  /I,  increases  with  beam  current  because  dcr/dC  is 

P  b 

2 

proportional  to  I^/a  ;  at  higher  values  of  1^,  the  effective  current  is 
frozen  in  at  an  earlier  time.  For  air,  and  most  other  simple  gases, 

1^  ~  lOkA  represents  a  transition  point;  higher  current  beams  are  mostly 
current  neutralized,  while  lower  current  beams  are  only  weakly  neutralized. 

If  1^  rises  to  its  full  value  over  a  time  long  compared  the  temporal 
delay  from  the  beam  head  to  the  pinch  point,  then  the  erosion  rate  and  In 
become  functions  of  the  beam  rise  rate  rather  than  the  peak  1^.  However, 
if  the  beam  propagates  far  enough,  the  rising  portion  of  I,  (t)  eventually 
erodes  away  and  Eqs.  (8)  and  (9)  become  directly  applicable. 


Gas  Conductivity 

We  have  already  had  a  number  of  occasions  to  refer  to  the  evolution  of 
a,  and  this  aspect  of  the  physics  also  has  a  major  effect  on  beam 
instabilities.  At  this  point,  we  shall  briefly  consider  the  principal 
mechanisms  that  underlie  plasma  conductivity. 

The  conductivity  may  be  written  neep,  where  u  is  the  electron 
mobility,  determined  in  general  by  both  electron-neutral  (e-n)  and 
electron-ion  (Spitzer)  collisions.  Frequently  a  dense  gas  is  only  weakly 
ionized  by  a  beam  and  e-n  collisions  dominate,  in  which  case  u  is 
independent  of  ng  and  only  weakly  dependent  on  temperature  Tg  (typically,  u 
«  1//Te>.  Thus,  a  depends  primarily  on  the  ionization  and  de-ionization 
processes  that  control  ng,  which  increases  rapidly  from  zero  to  ng  >  10 
In  addition  to  beam-collisional  ionization,  avalanche  ionization  (i.e., 
ionization  of  atoms  by  plasma  electrons  which  have  been  energized  by 
macroscopic  electric  fields)  may  occur,  particularly  near  the  pinch  point 


1 
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where  large  electric  fields  are  present.  De-ionization  is  usually  due 
primarily  to  recombination,  at  least  within  the  core  of  the  beam.  (In  some 
gases,  e.g.,  ,  attachment  can  be  important  at  large  radii  or  wherever  ng 

is  relatively  low.  This  will  be  ignored  here.)  We  may  thus  write  an 
equation  governing  the  evolution  of  a: 

dff/dC  =  a1(Te> +  a2(Te)a  -  3(Tg)a2.  (10) 

In  the  front  portion  of  the  beam  and  out  to  one  or  two  beam  radii, 

beam-collisional  ionization  [the  first  term  of  (10) J  usually  dominates, 

except  possibly  at  the  pinch  point.  In  this  case  (and  if  we  also  neglect 

the  ^-dependence)  the  radial  profile  of  a(r)  is  identical  to  that  of 

J^(r).  Moreover,  the  plasma  current  J^fr)  =  crE^  has  a  similar  profile, 

since  E is  only  weakly  dependent  on  r  within  the  beam.  [See  Eq.  (11) 

below.]  This  approximation  is  frequently  made  in  analytic  studies,  and 

greatly  simplifies  the  analysis.  Avalanche  [the  second  term  in  (10)]  has 

two  effects.  Avalanche  driven  by  E z  at  the  pinch  point  causes  the  central 

o(r)  profile  to  become  narrower  than  J^(r),  a  destabilizing  effect  that 

will  be  discussed  later.  Avalanche  driven  by  Ef  in  the  very  front  of  the 

beam  produces  a  very  broad  low-level  conductivity  out  to  the  large  radius 

b,  which  has  important  consequences,  e.g.,  the  large  inductive  logarithms 

in  Eq .  (8).  Further  back  in  the  beam,  where  recombination  [the  third  term 

of  (10)]  becomes  important,  <x(r)  becomes  broader  than  J^(r),  which  helps 

stability.  We  note  that,  in  a  complex  gas  such  as  air,  recombination  can 

depend  on  water  vapor  content  and  complicated  temperature-dependent 

12 

chemistry  effects,  e.g.,  the  formation  of  molecular  complexes.  On  a 
practical  level,  even  the  weak  dependence  of  u  on  T^  can  have  significant 
effects,  e.g.,  on  hose  instability,^  as  discussed  below.  To  model 
temperature-dependent  effects,  it  is  frequently  convenient  to  assume  that 
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(in  a  diatomic  gas)  Tg  is  determined  by  a  balance  between  ohmic  heating  and 
energy  loss  to  vibrational  modes  of  the  molecules.  In  this  case,  it  can  be 
shown  that  Tg  is  a  function  only  of  E/p,  where  E  is  the  electric  field  and 
p  is  the  gas  density.  We  can  then  write  the  coefficients  in  Eq .  (10)  as 
functions  of  E/p.  A  summary  of  these  coefficients  is  given  in  Refs.  13 
and  14. 

We  have  assumed  that  the  plasma  conductivity  is  local,  scalar,  and  is 
created  instantaneously  by  beam  collisions.  In  fact,  beam-gas  collisions 
do  create  high-energy  secondary  electrons  (delta  rays)  which  have  long  mean 
free  paths.  Although  relatively  few  in  number,  these  secondaries  can  have 
some  effect  in  spreading  and  delaying  conductivity,  and  in  responding 
nonlocally  to  electric  fields.  It  is  also  true  that  tensor  conductivity 
can  play  some  role  if  the  net  current  is  high  or  the  beam  radius  is  small, 
leading  to  strong  magnetic  fields. 


Beam  Instabilities 

As  a  magnetically  confined,  highly  ordered  system,  a  self-pinched  beam 

is  subject  to  a  number  of  instabilities,  which  are  driven  by  two  effects: 

(1)  If  there  is  a  substantial  return  current  I  s  -fl,  ,  then  there  is  a 

p  b 

repulsive  magnetic  force  between  1^  and  1^.  As  long  as  1^  and  1^  remain 
well  aligned  and  more  or  less  proportional  to  each  other,  this  merely 
weakens  the  pinch  in  an  orderly  way,  but  if  perturbations  lead  to  a 
separation  of  I  and  1^,  the  repulsive  force  can  drive  unstable  growth  of 
the  perturbations.  This  is  the  primary  mechanism  for  all  of  the 
instability  modes  except  hose.  Each  has  a  threshold  value  of  f  below  which 
the  mode  is  stable. ^  (The  hose  mode  is  unstable  even  if  f  =  0,  but  is 
further  destabilized  by  nonzero  return  current.)  (2)  Even  if  I  =  0, 
symmetry-breaking  distortions  can  drive  locally  destabilizing  magnetic 


forces.  Finite  plasma  resistivity  plays  a  key  iole  tieie.  Because  a  t  °>, 
magnetic  field  lines  are  not  frozen  into  the  plasma,  and  instabilities  can 
occur  on  the  time  scale  for  beam  motion,  rathei  than  the  much  slower 
hydrodynamic  time  scale  for  the  plasma.  But  because  o  *  n,  the  field  lines 
are  subject  to  destabilizing  phase  lags  as  they  tiy  to  follow  beam 
distortions.  This  mechanism  particularly  drives  the  resistive  hose 
instability,  which  is  the  most  notorious  of  the  beam  instabilities. 

The  linearized  normal  modes  of  a  beam  can  lie  chara<  tetired  by  a  pait 
of  quantum  numbers  (m,n),  where  m  indicates  @- dependence  exp(imA),  and  n  is 
the  radial  mode  number,  roughly  speaking  the  tuimbet  of  oscillations  within 
the  beam  radius.  The  first  few  modes,  shown  in  Fig.  1.  ate  the  most 
important  for  pinched  beams:  (m  -  0.  n  -  1)  is  the  sausage  »f"d<-,  toughly  a 
self-similar  expansion/contraction  of  the  beam:  (m  -  ",  u  .' )  is  the 
axisymmetric  hollowing  mode,  in  which  the  beam  density  alternately  hollows 
out  and  peaks  on  axis;  (m  =  1,  n  =  0)  is  the  hose  mode,  in  which  the  beam 
thrashes  around  more  or  less  like  a  firehose,  without  a  great  deal  of 
internal  distortion.  The  f ilamentation  modes  (m  >  2  oi  n  >  1)  generally 
have  a  high  threshold  value  of  f^  and  are  believed  to  he  stable  for 
pinched  beams,  except  in  annular  (usually  totaling)  beam  equ i 1 i hi i a . 

Axisymmetric  Hollowing  Instab i  1  i  ty 

The  axisymmetric  hollowing  instability  was  dismveted  in  <omputei 
simulations^  as  a  particularly  violent  instability,  leading  to  tapidly 
growing  radial  oscillations  that  destroy  the  beam  only  a  few  nan.ise. ends 
behind  the  pinch  point.  Computer  simulations  have  piovided  a  detailed, 
quantitative,  and  rather  surprising  picture  of  its  natmo.  Figute  shows 
the  growth  of  the  instability  as  a  function  of  z  and  C  (used  as  independent 
variables  in  place  of  the  usual  z  and  t  ),  and  Fig.  )  shows  the  radial 
profile  .  I  b ( t  ,  C)  characteristic  of  the  instability. 


Since  the  simulations  show  the  instability  in  the  large-amplitude 
nonlinear  stage,  where  large  radial  oscillations  are  apparent,  it  was 


k.v 


thought  at  first  that  this  was  basically  a  sausage  instability,  although 
hollowing  of  the  beam  profile  is  apparent  in  Fig.  3.  However,  when  the 
beam  profile  was  constrained  to  a  self-similar  shape,  the  instability 


disappeared,  thus  indicating  that  hollowing  is  an  essential  feature 

16 


Furthermore,  a  linearized  analytic  theory  of  the  sausage  mode  shows  that 
instability  is  not  expected  for  a  beam  injected  into  neutral  gas. 

[Basically  this  is  because  beam  perturbations  create  similar  conductivity 
perturbations  through  Eq .  (10).  This  effect  inhibits  spatial  separation  of 
1^  and  I  ,  and  thus  reduces  the  growth  rate  for  all  modes;  it  just  barely 


suffices  to  stabilize  the  sausage  mode.]  Furthermore,  when  the  avalanche 
term  in  Eq.  (10)  was  artificially  turned  off,  the  instability  disappeared. 
Finally,  it  was  found  that  when  E/ p  was  extremely  large  (specified  below), 
the  instability  turned  off. 

The  essence  of  the  hollowing  instability  is  as  follows.  Behind  the 
pinch  point,  where  Ampere's  law  (in  its  axisymmetric  form)  is  valid,  one 
can  show  that  the  electric  field  has  the  radial  profile 


E  (r)  «  1 n [ ( 1  +  b2/a2)/(l  +  r2/a2)]. 


(11) 


Ez  has  a  weak  (logarithmic)  maximum  on  axis,  which  is  usually  ignored  in 


analyses.  However,  over  a  wide  range  of  parameters,  the  exponentiation 
rate  for  avalanche  ionization  rates  increase  as  a  high  power  (4  to  6)  of 


E/p.  Thus,  ct,  and  also  the  plasma  return  current  J  =  E  o,  become  strongly 

P  z 


peaked  there,  and  the  repulsive  magnetic  force  hollows  out  the  beam  and 
blows  it  out  to  large  radius.  This  expansion  reverses  itself  only  because 
the  radial  expansion  of  the  beam  current  decreases  the  inductance  L  of  the 
system.  Roughly  speaking,  LI  tends  to  be  constant,  so  the  defocusing 
plasma  current  decreases,  or  may  even  reverse  itself  so  as  to  augment  the 


pinch.  ’  The  beam  current  then  comes  crashing  back  onto  the  axis,  and 

the  cycle  is  repeated  with  rapidly  increasing  amplitude. 

The  key  quantitative  features  revealed  by  the  simulations  are  that  the 
instability  occurs,  over  a  wide  range  of  densities,  only  if  two  conditions 
are  met.  First,  f  >  50%,  in  qualitative  agreement  with  prior  analytic 
predictions.^  Ibis  is  easily  understood:  even  if  the  return  current 
flows  in  a  profile  that  is  very  narrow  compared  to  the  beam  current,  the 
pinch  is  destroyed  overall  if  and  only  if  >  J^/2.  Secondly,  the  value 
of  Ez/p  at  the  pinch  point  (in  air)  must  fall  into  the  range 

13  MV/m-torr  <  E^/p  <  50  MV/m-torr.  (12) 

The  lower  limit  on  E^/p  ensures  that  avalanche  is  strong  enough  to  play  a 
significant  role.  The  upper  limit  is  due  to  the  fact  that,  although 
avalanche  is  very  strong  at  large  values  of  E/p,  it  no  longer  increases 
rapidly  as  a  function  of  E/p. 

In  order  to  avoid  the  hollowing  instability,  it  is  thus  necessary  to 
keep  the  maximum  value  of  E z  below  the  lower  limit  of  Eq.  (12).  (In  gas  at 
any  significant  fraction  of  atmospheric  density,  the  upper  limit  is  not 
exceeded.)  This  can  be  accomplished  in  several  ways:  (1)  By  limiting  the 
rise  rate  of  I^(t),  so  that  the  current  at  the  pinch  point  is  not  above  a 
critical  value.  (2)  By  limiting  the  peak  value  of  1^.  Even  if  I ^ ( t )  rises 
instantaneously,  instability  occurs  (in  air)  only  if  1^  >  20kA  multiplied 
by  the  density  in  atmospheres,  for  beam  radii  in  the  vicinity  of  1cm.  (3) 
Increasing  the  beam  radius. 

These  conclusions  were  subsequently  tested  in  an  experiment  performed 
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on  the  IBEX  electron  beam  facility  at  Sandia  National  Laboratories.  The 
hollowing  instability  was  clearly  seen  to  occur  at  air  densities  below  80 
torr,  in  quantitative  agreement  with  predictions  based  on  Eq .  (12),  and  to 
turn  off  at  higher  pressures  where  E  /p  became  too  small. 


Hose  Instability 

The  resistive  hose  instability  is  the  most  important  impediment  to 
propagation  of  pinched  beams.  It  is  observed  in  nearly  all  beam 
propagation  experiments,  and  has  been  studied  extensively  by  means  of 
linearized  analytic  theory  as  well  as  both  linearized  and  nonlinear 
numerical  simulations.  The  analysis  indicates  that  the  instability  is 
always  present  for  a  pinched  beam  injected  into  neutral  gas,  but  that  it 
can  be  minimized  by  limiting  the  beam  duration  and  reducing  the  level  of 
initial  perturbation. 

Analytic  studies  of  the  hose  instability  have  usually  been  based  on 
linearized  theory  for  perturbations  to  an  axially  uniform  beam  equilibrium, 
i.e.,  J^o(r,0  independent  of  £.  Hose  normal  modes  take  the  form 
f(r,C)exp[i(9  +  Qz/c)].  It  may  also  be  assumed  that  the  equilibrium  plasma 
current  JpQ(r,C)  and  (in  the  simpler  theories)  conductivity  aQ( r,C)  are 
C-independent ,  in  which  case  the  normal  mode  dependence  reduces  to 
f(r)exp[i(0c  +  wC  +  Sz)/c]  and  a  dispersion  relation  w(S2)  is  sought. 

Effects  associated  with  the  beam  head  and  pinch  point,  e.g.,  low 
conductivity,  space  charge,  and  incomplete  pinch,  are  usually  neglected, 
and  in  this  spirit  the  electrodynamics  are  calculated  simply  from  Ampere's 
law  (exceptions  are  Refs.  20,21). 

The  earliest  version  of  hose  theory  (the  "rigid  beam"  model)  assumed 
in  addition  that  the  perturbation  of  each  "slice"  of  the  beam  consists  of  a 
transverse  displacement  by  an  amount  Yexp[i(«C  +  S?z ) J ,  with  no  internal 
distortion.  As  a  nearly  exact  consequence  of  the  linearized  Ampere's  law 
the  vector  potential  A  (r,z,C)  is  also  displaced  from  the  axis  of  symmetry, 
by  an  amount  Dexp[i(wC  +  S5z)|,  without  internal  distortion.  The  problem 


(13) 


2  2  2  2 
a  Y/3z  =  52-  (D~Y)/c 

P 


d  +  cx^aD/ao  =  y, 


2  2 

where  x^  =  itaa  /2c  is  the  "dipole"  magnetic  decay  time.  Equations  (13) 


and  (14)  lead  to  a  dispersion  relation 


2  2  2 

ioox,  =  -52  /(k-  -Q  ), 
■L  P 


which  correctly  shows  that  oscillation  in  z  scales  to  A„  while  growth  in 

P 


C  scales  to  the  x,  .  However,  the  infinite  growth  rate  predicted  at  52  -»  k_, 
i  P 


indicative  of  absolute  instability  in  the  beam  frame,  is  incorrect. 


The  crucial  oversimplification  in  Eqs.  (13)-(15)  is  the  implicit 


assumption  that  all  beam  electrons  oscillate  at  a  single  resonant  betatron 


frequency,  i.e.,  that  the  potential  well  pinching  the  beam  is  simply 


harmonic.  This  would  be  true  for  a  flat-topped  current  profile,  but  the 


rounded  profiles  of  and  introduce  anharmonici ty ,  and  therefore  a 
dependence  of  k^  on  the  amplitude  of  an  electron's  orbit.  When  this 


feature  is  introduced  into  the  modeling,  as  was  first  done  by  using  the 


"spread  mass"  formalism,  the  dispersion  relation  exhibits  a  finite 


maximum  growth  rate.  For  example,  Lee  finds 


-iwXj  =  3x^-6x^  +  6(x^-x^)[in  +  ln(l/x^-l)], 


2  2 
/  r\  *- 


where  x  s  52  /  52^  .  Equation  (16)  is  illustrated  in  Fig.  4,  for  a  beam  wii 


no  return  current.  Most  significantly,  the  instability  is  convective 


backward  in  the  beam  frame;  hence,  it  reaches  a  maximum  amplitude  at  any 


given  point  in  the  beam  and  then  decays.  Thus,  the  hose  amplitude  can  be 


limited  by  limiting  the  beam  duration  to  a  few  growth  lengths  and  by 


insuring  that  the  beam  is  initially  quiescent,  so  that  hose  modes  have  to 


e-fold  many  times. 


s 


Subsequent  to  Lee's  pioneering  work,  hose  modeling  has  been  extended 
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in  many  ways.  More  sophisticated  macroscopic  models,  linearized  Vlasov 

it.-  24  .  .  ,  17,18,23-25  .  .  .  _  ,  „  , 

calculations,  and  simulations  have  been  used  to  treat  beam 

dynamics.  Plasma  return  current  has  been  included  and  found  to  be  strongly 
23  24  26 

destabilizing.  ’  ’  Self-consistent  treatment  of  conductivity  evolution 

introduces  a  variety  of  effects.  A  fully  analytic  linear  theory  has  been 

developed  which  includes  the  linear  increase  of  equilibrium  conductivity 

ct  ( r,C)  with  C  due  to  beam-collisional  ionization  of  the  gas;  in  this 
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case,  instability  grows  as  a  power  p  of  C»  rather  than  exponentially. 

Since  p  is  inversely  proportional  to  a  ,  growth  is  most  rapid  just  behind 

rhe  pinch  point.  (The  formalism  is  invalid  ahead  of  the  pinch  point). 

Furthermore,  dipole  perturbations  of  the  beam  induce  dipole  perturbations 

of  a  through  Eq.  (10).  This  significantly  reduces  the  growth  rate  in  the 

presence  of  plasma  return  current,  particularly  for  low  frequency  modes,  by 

2  A 

inhibiting  separation  of  the  beam  current  from  the  return  current.  One 
consequence  of  this  is  that  the  hose  growth  rate  typically  decreases  with 
increasing  1^,  as  shown  in  Fig.  5.  This  is  the  net  result  of  three  pieces 
of  physics:  (i)  a  1^,  which  favors  higher  currents;  (ii)  the 
destabilizing  effect  of  current  neutralization,  which  favors  lower 
currents;  (iii)  the  effect  of  dipole  conductivity,  which  swings  the  balance 
to  higher  currents.  Even  the  T^-dependence  of  a,  i.e.,  the  inverse 

13 

dependence  of  a  on  E/p,  can  be  included  in  a  fully  analytic  theory;  since 
Ez  is  largest  at  the  pinch  point  and  steadily  decreases  thereafter,  this 
further  accentuates  the  tendency  for  hose  to  grow  rapidly  at  the  pinch 
point  and  very  slowly  further  back  in  the  beam. 

Numerical  simulation  has  been  essential  to  detailed  understanding  of 
hose.  It  permits  self-consistent  treatment  of  the  radial  and  C 


dependence  of  beam  equilibria,  of  phenomena  near  the  pinch  point  where 


O.AV.V.V.V.V. 


Ampere's  law  is  invalid,  and  of  nonlinear  effects.  Conversely,  the  hose 

instability  has  stimulated  the  development  of  innovative  simulation 
24  25 

models.  ’  One  conclusion  from  these  studies  is  that  avalanche 
ionization  at  the  pinch  point  is  lethal,  strongly  driving  hose  as  well  as 
axisymmetric  hollowing,  and  for  similar  reasons. 

Conclusions 

We  have  seen  that  propagation  of  beams  in  dense  gases  is  limited  in 
varying  ways  by  energy  loss,  radial  expansion,  nose  erosion,  and  hose  and 
hollowing  instabilities.  For  the  most  part,  these  limitations  are 
minimized  by  going  to  higher  energies  and  currents,  fatter  beams,  shorter 
pulses,  and  beams  which  are  more  quiescent  at  injection. 
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Beam  median  radius  R, .«  as  a  function  of  C  and  z.  Case  (a)  is 
stable,  and  shows  only  the  expanded  and  slowly  eroding  beam  head. 
Case  (b)  is  unstable,  and  shows  the  large  radial  oscillations 
characteristic  of  the  hollowing  instability.  (From  Ref.  14.) 
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Fig.  3.  Beam  current  density  as  a  function  of  C  (the  distance  behind  the 
beam  head)  for  fixed  z,  showing  hollowing  and  peaking.  (From 
Ref.  14.) 
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Peak  hose  growth  rate  as  a  function  of  beam  current.  The  curve  is 
an  analytic  result  from  Ref.  24,  while  the  dots  are  the  result  of 
simulations  with  the  mul t i - componen t  code  VIPER,  also  taken  from 
Ref.  24. 
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